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ABSTRACT 


Demand during lead time is considered as the sum of a random number 
of random variables. Discrete random variables are considered as ds- 
scribing demand and lead time and probability generating functions are 
used to determine the probability generating function of demand during 
lead time, The probability distribution of demand during lead time is 
extracted from the probability generating function by using a power 
series expansion and the calculation of probabilities from the expansion 
are discussed, The mean and variance of demand during lead time are 


also obtained from the probability generating function. 
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Ly Introduction 

The two major problems in any inventory control system are to 
decide upon the amount to order each time replenishment is necessary 
and to decide what should be the minimum stock level, of any item, 
at which an order should be placed, The second problem is the subject 
of this paper. 

The minimum stock level at which an order for replenishment of an 
item is the reorder point, denoted by R. The reorder point is a funce 
tion of the demand and the lead time of an item. The demand for an 
item is the rate at which the stock of the item is depleted, The lead 
time of an item is the length of the interval from the time a replenish- 
ment order is placed until the order is in stock, ready to satisfy the 
customer's demand, 

A simple deterministic model of an inventory system with constant 


demand and lead time is shown in figure l. 
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Deterministic Inventory Model 


Figure 1 








The stock is depleted at a uniform rate from a stock level of S 
until it reaches the reorder point, When the stock level of R units 
is reached, a replenishment order for S units is released and the 
order is received at time A, just as the stock is exhausted, The 
stock replenishment raises the stock level back to S units and the 
process is repeated. 

Rarely is demand uniform, however, and there is uncertainty about 
when the R units of stock on hand at the time of placing an order will. 
be exhausted. A more realistic model of an inventory system with 


variable demand and constant lead time is shown in figure Ze 
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Variable Demand Inventory Model 


Figure 2 


S, R and A are the same in figure 2 as they were in figure l. 
The expected depletion of stock is along the line S A. If the lead 
time was deterministic, an order placed at R, time X, would be rsq= 
Gaived at time A. In the time interval during which replacement 


stock is an order, K A, there can be a variation in the demand; 





therefore, the demand may deplete the R units of stock along any line 
between a maximum demand, R A', and a minimum demand, R A", shown by 

the shaded area A! 2 A", The variation in demand allows for ths pos- 
sibility of depleting the stock at time A' and not being able to meet 
demands for stock in the period A' A, There is also the possibility 

that the demand could be as low as R A" and in this case there would 

be B units of stock on hand when the replacement order is received, 

A third and the most realistic model, considering both variable 
demand and lead time, is shown in figure 3. The first two models 
were special cases of this model and therefore model 3 will serve to 
describe the inventory system considered throughout the rest of this 


paper. 
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Variable Demand and Lead Time Model 


Figure 3 


The variation in the lead time will cause the replenishment 


time interval X Y to vary, Time X is the time at which the reorder 





point, R, is reached, but it is no longer certain that the order 

will be received at time Y. The order may be received at time U or 
time V or any time in between, if we let time U and time V represent 
the minimum and maximum times of receipt of the order, Figure 3 is 
similar to figure 2, but the zero stock level of figure 2 is ths 

same number of units below R as G is in figure 3. Figure 3 then is 
the same diagram except that the entire stock is increased by G units 
to allow for the variation in lead time as well as the variation in 
demand. Let us consider the four possible extreme cases of maximum 
and minimum demand and lead time, A replenishment order is placed 
when the stock on hand is R units at time X. If the maximum demand 
and the minimum lead time occur, the stock will be depleted along 

line R A until there are D units at time U when the replacement stock 
is received, If the minimum demand and the maximum lead time occur, 
the stock will be depleted along line R A" until there are E units in 
stock at time V when the order is received, Minimum demand and mini-= 
mum lead time will result in F units of stock on hand at time U when 
the order is received, In each of these three cases, R units of stock 
have been more than enough to provide protection against stockouts 
with both variable demand and lead time, Only in the case of maximum 
demand and minimum lead time would the stock be below © units, The 
last case is the combination of maximum demand and maximum lead time, 
The stock will be depleted along the maximum demand line R A‘ and will 
be exhausted at time Z. The order will not be received until the maxi- 
mum lead time, time V, therefore the system will be out of stock for a 


period of V-Z time units. 





From figure 3, it can be seen that if the demand follows the 
line R A and time Y is the lead time, then the demand during lead 
time is equal to R = C units of stock, It can also be seen in 
figure 3 that as the height 0 C is increased the intersecticn of 
the actual demand and the actual lead time, whatever they may be, 
will intersect above the zero stock level more often than below, 
The reorder point, R, is then the sum of the expected demand during 
lead time plus a safety level or buffer stock of C units, The re=- 


order point equation is then; 


Ren Set (1.1) 
where 

R = reorder point 

Z= expected demand during lead time 


VS L= variable safety level 


The safety level is denoted as a variable safety level because the 
protection against a stockout increases as the number © is increased, 
The safety level can be varied so that no stockouts occur or soms 
small number of stockouts occur in a specified number of replenish-= 
ment periods, That is, the inventory manager may decide that a 
safety stock large enough to protect against no stockouts would be 
too expensive and therefore he is willing only to maintain a safety 
stock that will protect against stockouts ninety or ninety-five per-= 
cent of the time, The safety level is also a function of the reorder 
quantity since the larger the reorder quantity for a given lead time 


the less frequently the system is exposed to a stockout. This 





consideration is taken into account when the protection level is 


determined and does not affect the reorder point equation GE alle 





2. Discussion of the Problem 

The availability of alternative courses of action implies the 
need for managers to make decisions, The selection of the reorder 
point for an item in the manager's stock ig then a desision problsm, 
The reorder point rule must satisfy an objective function (such as 
maximum service) subject to some constraints imposed by higher 
authority, available funds, available facilities, etc, In almost 
all cases, the method by which the problem is solved requires a model 
of the operations. Such a model is based on a set of relations link-= 
ing the important variables. Figure 3 is the geometric modséi or ths 
inventory system being considered and the reorder point equation (1.1) 
will lead to the mathematical model which will be used to determins 
the probability distribution function of demand during lead tims, 

The reorder point equation is a rule by which the inventory 
manager may decide when to reorder replacement stock, The variable 
Z is the mean of the probability distribution of the demand during 
lead time and the variable V S L is a function of the variance of the 
same probability distribution, Therefore; in order to use the reorder 
point equation, the probability distribution function of demand during 
lead time or at least the mean and variance of the distribution func= 
tion must be known. 

Let Z be a random variable which describes the demand during 
lead time, then the number of demands that occur during a lead time 


may be written as: 


Z = D. = Do + eo90 she (2.1) 





where 
Z = the number of items demanded during a lead time 
D,= demands : 1 = |gmeos 5 


L = lead time 


It is assumed that the demands are independent and identically dis+ 
tributed and demand and lead time are independent discrete random 
variables with known probability mass functions. In general, this 
paper will only consider discrete random variables sinca demands 
naturally occur at discrete intervals and lead times are normally 
described in days, weeks or months, Equation (2.1) is then a model 
which describes the demand during lead time by linking together ths 
random variables D and L and is a model that can be used to detere 
mine the expected demand during lead time, vA. and the variable 
safety level, VS L. Equation (2.1) can be described more concise 


ly as: 


L 
= >, Gere 
Z = 4 D, (2.2) 
Using equation (2.2) as the model of demand during lead tims, 
it, will be shown that the mean and variance can be computed for any 
probability mass functions that describe demand and lead time or 
from the mean and variance of demand and lead time if that is all 
that is known about them. It will also be shown that the distri- 
bution of demand during lead time can be found, in the discrete 
case, by the use of probability generating functions, After the 


probability distribution is found, it will then be shown that it 





can be used to set the reorder point using a variable safety level; 
given the level of customer service desired, The final section 
treats a special case in which lead time is a continuous random 
variable and the resulting distribution of demand during lead time 


has a familiar negative binomial form, 





3, The Mean and Variance of Demand During Lead Time 

The mean and variance of Z are derived using general probability 
mass functions for D and L. 

Let: 

p (a) be -the mass function of demand with mean D and variance 


Oo? 


D 


p (1) be the mass function of lead time with mean L and vari- 
2 
ance OC, 


p (z) be the mass function of demand during lead time with 


= 2 
mean Z and variance OC, 


The mean of Z is the expected value of Z and is denoted as Z 


if 


E [B(2/2)| 
g [L £(D) 
E(D) E(L) 


E(Z) 


= bee 


Therefore the expected demand during lead time is 


fc 18) 1 


The variance of demand during lead time is derived as follows: 


Var (Z) 


E [var (2| 1) + Var [ B(211)| 


E E — + Var E BD) 


Var (Z) 


10 





= E(L) on + poy 0, 


¥ 


The mean and variance of demand during lead time are functions 
of the means and variances of the demand and lead time only, there- 
fore a knowledge of the mass functions of demand and lead time is 
sufficient to calculate the mean and variance of demand during lead 
time, They may also be calculated empirically using historical de- 
mand and lead time data to estimate the means and variances of demand 
and lead time, The mean and variance of a distribution function may 
also be calculated from the probability generating function of the 
distribution, The probability generating function of a probability 
distribution function is discussed in the next section and the calcu 
lation of the mean and variance from a probability generating func= 


tion is shown in Appendix I. 


aa 








4. Distribution of Demand During Lead Tims 

In this section, the probability generating function of demand 
during lead time will be derived in general. The generating funce- 
tion will then be used to obtain the probability distribution of 
lead time demand by expanding the generating function in a power 
series and extracting the coefficients, ‘frobability generating 
functions are defined for discrete random variables only. There-~ 
fore; the discussion will be limited to discrete random variables, 

The probability generating function of demand during lead time 
is the generating function of the sum of a random number of random 


variables and is defined as 
OO 
= nie = 2y st > [ed 
h(s) =E(s ) = =i Berge“ 


To derive h(s), let the probability generating function of demand 


and lead time be 


e(s) = E(s”) 
f(s) = E(s>) 
and 
h(s) = E(s2) = E (3/1 = 1) 
now 


Z D ! 
E(s IL =1) =a 5 woof) 
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cp 
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since the Dj; were assumed independent and identically distributed, 


Also since f(s) = B(s?), we have 


E(s“|L = ee [e(s)]* 


and 
B(s“|L) = [#(s)} : 
E |e(s4 |x) | = B 4 (els) ul 
= g le(s) 
since 
g(s) = E(s”) 
»°. h(s) =g¢ (f(s) 


Thus, if the distribution functions of demand and lead time are known 
or assumed, the probability generating functions of demand during lead 
time can be derived, Then with a knowledge of h(s), the probability 
distribution of demand during lead time is known. Examples showing 
the computation of the distribution of Z will be given later, 
Implementation of the model given in equation (2.2) may generally 
take place using familiar discrete distributions for the variables, 
Since such distributions have convenient mathematical properties. 
These familiar distributions usually have positive probability at zero 
and this situafion usually is not compatible with lead time. Prac= 
tically, lead “time can never be zero, if for no other reason than it 


takes some time to prepare an order, On the other hand, if the data 
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in the inventory control system is based on the month as the basic 
time unit, a lead time of one day could be considered zero. In any 
event, if the mean of the lead tims distribution is small, then ‘ths 
use of a popular distribution will rasult in probability at zero which 
is not negligible. Mathematically, the model is correct since when L 
equals zero the sum is empty and Z is zero, Demand during lead time 
can also be zero when demand is zero during the lead time or when 

both demand and lead time are zero. 

There are several methods by which the input distribution can be 
adapted to consider either allowing or not allowing some probability 
for lead time to be zero, In the first casa, when lead time can be 
zero, the straight forward use of the model with distribution fune- 
tions of demand and lead time that arse defined on the non-negative 
integers will have provided the desired results, There ars two methods 
that may be used to prevent the consideration of zero lead time in the 
model, These methods redistribute the probability mass functions from 
a range over the non-negative integers to a range over the positive 
integers. The first is translation, i.¢., let L = 1+X, where X takes 
non-negative integral values; and the second is truncation, i 4., 
use a frequency function that is conditioned on ths variable being 
positive. Let us consider the Poisson distribution as an example, 

The Poisson distribution is usually defined as: 


ere es 





p(x3A) = LOL Og Ore oo ec 


The mass function can be translated over the positive intsgers by 








writing it as: 


an x1 


paA) = ie 


, for A> 0a Cae, 


or the distribution function can be truncated by using a conditional 


probability, That is; 











=A x 
p [x = x|x>0] =2 A 
x!(tee ~) 
x 
~~ D 
e = x! 
: ee 
on) SS tor Ae Ox = ieee 
(e%-1) x! 


There are then three choices offered for adjusting the input to fit 
the situation or the preference of the user. 

Let us now consider some generating functions of demand and lead 
time so that they will be available for use in determining the prob- 
ability generating function of Z. The Poisson and negative binomial 
distributions have frequently served as good approximations for the 
distribution of demand so they will be used in the following examples. 
lead time distributions have not been as well documented, but the 
negative binomial distribution is a good model for discrete waiting 
time random variables and when the parameter r equals one, the nega-= 


tive binomial becomes the geometric distribution. As a second lead 


ike. 








>_> 
= = ==) 





time distribution, the Poisson distribution will be used. The use of 
the Poisson distribution for lead time may not ba a good sxample, but 
with a large mean, the Poisson distribution roughly approximates a 
normal distribution with the mean equal to the variance [a]. 

The generating functions for the Poisson and negative binomial 
distributions in sach of the three cases are given below, 
For the Poisson, when 


wae 
2c 





p(xsA) = 


The probability generating function is: 


CO 
fle) = (3) = Y. 3* Pix = 
feo, 
= ee) (a4) 
when » P 
— Xo 
p(xs A) = 2A » X= 1, Pate ae) 
A(s=1 
HOR eS (4.2) 
when 
d ) A* 
p(x ae ea to. 
9 (e-1) 7 9 9 bY) 
\s , 
d 
e <= 


for negative binomial distribution, when 
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+x=1 
pi snap a= (* . p’ q* Sc ON ipl 
r 
£(s) = ( : (4.3) 
i-q s 


when the distribution is translated 


+X=m = 
p(xsr,p) = (" ‘i _ pr gX-1 Fax. 1,02 ae 


bidiea, re ) 


when the distribution is truncated 


ee os 
x p q 


q=p 





p(xsr,p) = 9X = 1, 25 co 


As an example of determining the probability generating func- 


tion of demand during lead time, let 


D be a Poisson random variable 
L be a translated negative binomial 
then 


f(s) =. A(s-1) 


and 


2 
o~ 
63) 
—— 
il 


g f (s) 


sk 











Om 


= pt he e Ms-1)k ee og k-1 


k=1 k=1 


9:8) 
r _ A(s-1) a A(s-i) (k=1) frtice2\ — k-1 
p 8 k=4 © = q 


a0 

ayy kat 
r , A(s-1) >. frtt=2) A(s=1)} 
RS k=1 k=1 


of g A(8-1) [ye Ale=1) |™ 


The probability generating functions of demand during lead time, Z, 
were derived since they can be used to find the prcbability of de= 


mand during lead time. That is, since 


OD 
h(s) = E(s¢) = >, s* Pp iz ~ | 


Then in the power series expansion of h(s), the coefficient of s* 


is the plz = ¢). For example, writing out several terms of h(s), 


we have: 


h(s) = Piz = o]+ s Pk = 1] +o? Plz = 2] +... 


It therefore is necessary to put h(s) in the form of a power series 
in order to determine the probability of demand during lead time 
for any k. 


Using the probability generating functions for the variables Z 
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in the previous examples, we can put them in a more useful form. 


For Poisson demand and translated negative binomial lead time 


h(s) =e N(8-1) mere | 
1-9e A(s-1) 


(oe) 


» A(s-1) or >. (rei) fag a(sot)]" 


k=0 


QO 
omer a ee k os A(kt+H1) - (k+1) 
ss 2, ( 2 


it 
Ts 
M1 


oo 
Gow gk 9 A(k+1) >. 5f Later 


n=0 


ey git en A(k+1) LA(k+1)] ® 
n! 
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é me [2 = a| is the coefficient of s" in the power series 


and 


Bey A(k+1) [A(k+1)]" 


“eee” 
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k=0 n! 
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The calculation of these probabilities can be reduced if the sum 


from zero to infinity is reduced to a finite number of terms with a 


very small remainder, A method for finding the number of terms to 
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be summed until the remainder is negligible will be shown along with 
a computing technique for determining the probabilities, 


To find a bound for the remainder of the series, let 


K 


2. (eve (qe~* )* (k+1)" 


Av = biqe = 
RK” KOK Ke0l k 


co 
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and 
CO Go 
a — S_OK HL 
Ceca OT imo OR 
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K+] = rttet) gktttl roa ie 
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= (r+K+1)1(K+1)! (a a! 


(K+1+1)1(r+K)! K+2 


— _(r+K+1) (r+K+1-1)... (K+l+r-(r-2) ) (28)" ol 
(Ktr) (K+r-1)... (K+tr-(r-2) ) K+2 


_ (r+kt+1) (r+K+1-1)... (K+1+2) (222)" ral 
(Kt+r)(K+r-1) ... (K+2) K+2 


Dividing term by term 
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It is possible to dominate the right hand side with a integral pro- 
vided the integrand is monotone decreasing. This will be true if K 


is sufficiently large, Thus, let 


-B=lna=-=-A+i1nq 


and 
1 nt+1 7 
¥ (2) = (1 +e o7B1 
K 
then 
1\"2 _ a n+1 
va) = (+=) fa [- (1+—) 5+ 
K K K 


Setting the derivative equal to zero, we have 


. » the maximum for 1 eccurs if K > a 


hencs we must use K such that it is at least this large. 


Now for Kk > Se 
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pa) ac 

é: 1 \ +1 1 >» 1 \ n+ 1 
2, (1+—] a = 1+ (1 + — a 
1=0 K 


making a change of variables, let v=Bx, then dv=Bdx, and 


GO e n+1 oo 
a n+1 1 
(+B 2 2) te | dere 
1=0 j=0 K/ 8B 0 
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which is a bounded function of K, thus 


R, S Uri [ + Amt ( yn] 
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now if KB >n+1 


22 











v frtkK x + 
Ry < a eae A xaey | } tants ° | 


In practice, ail of the constants and parameters in the above 
equation are known except K. K can be determined by setting the 
equation equal to a very small number and solving for K. The 


error in P (Z=n) will then be less than that small number times 


sks re 
ie’ ac, 
nl 


Now that the value of K has been approximated, the probability 


distribution of demand during lead time can be found as 
K 


n 
p[z=n}=e" zt 


based 
ni k=0 K . 


= (+1) 
In order to simplify the calculation of this number manually or in 


a computer, the technique of "nesting" will be shown, 


K 
S. = aX(—44)% = q4{T) go (tT) 22 394 208 4 ee aS (K+1) 
a4 1 2 K 


K 
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= ( + ra 29 (+ ae - 7 (1+. -.(14 — eal ))--}) 


This can then be written in the compact form, 
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Let 
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and 
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B. = V+ A 

J j j= 
where 

+K=j—1 K=-j+1 : 

4 (— and B= 1 

J Kaj K-j : 
then 
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As another example 
Example 2 
let 
D be a negative binomial random variable 


L be a translated Poisson random variable 





Then 
r 
f(s) -(+ from eq. (4.3) 
i-qs : 

and 

e(eja=.s epee!) from eq. (4.2) 
Then 

h(a = 2 (a) "¢ | e(s)] = f(s) orle(s)-1] 


al, 


Tet us again find the Maclaurin expansion 
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To find a bound for the remainder of the series 
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We then have to choose minimum K such that ¥<1 so the resulting 


geometric series can be summed, i.e., 


Ye i( ° 2) | = 
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and to solve for K we have 


ak gf a b ( 1 ) = [x (K+1)+n-1] ! prk e7Agh py 
“Me ¥ [r (k+1) -1)int [- pr (1+ op FJ 


Upon setting this equation equal to an allowable srror, K can be 


determined , 


The estimate of the remainder is the value of the allowable 
error from which K was determined, 


Then, given K 


A Sa were (\ pr)* 
K 0 n k! 
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The “nesting” technique is helpful in computing the probabil- 


ities of demand during lead time, 


p [z= n| = en4 pf a” ay , T 21 
and 
Ay = (r+n-1)! a (2r+n-1)!(Ap_) 4 (3r+n-1) 1 (Ap™)? aes. 
(r=4) In! (2r-1) In! (3r-1)! n! 
. (rtn-1)! ( . (2a emer) (dp) (1 ie ene 
(r=1) In! (2r-1) 4... (r) 
“7 am (rk+n—1) 000 (rk+n-r) (A pt) (1 As Borns 
(rk=1) ... (rker)k 
saa , _(rktn-1) ooo (rk +n-r) (A p’ ) ))++)) 
(rK=1) e@oo (rK-r ) 
Thus, if 
—_ ee oe ee ooo sagt pare r 
3 =| ? =e I “i Tey | 5 ‘ r(K-j)-r AP 
and 
Co = a 9 then 
Cai = 1+ B, C, 
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and 


AK = eeeiirie ( + By-2 Cx-2 | 

Now that the distribution of demand during lead time and methods 
for calculating the probabilities have been discussed, it is a simple 
matter to find the number of items of stock that are required as the 
reorder point, The reorder point must be the number of units of stock 
that will last until new replacement stocks arrive with the probabil- 


ity of a stockout during the lead time small. 
Let a = Acceptable probability of a stockout during lead time 


Then 
Pp [z>z]<e 


must be solved for the value of R, 


1-p[z>R| = P[z =o] +P [z= 1] +...+P([z=R| 


e o Ris the smallest number such that 


ea 


n=0 
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5. Negative Binomial Distribution of Demand During Lead Time 


It would be very helpful if the distribution function of de- 


mand during lead time could be found for any distribution functions 


of demand and lead time in a familiar form, but this is not possible. 


It can be done in some special cases. One known case is given belo 


as an example, Consider deriving the mass function of demand during 


lead time where: 


(1) 


(2) 


The demand is random and the number of demands d 
during a fixed interval of time t has a Poisson 
distribution, 


-at 4\d 
a Ed. 


gd = 051525006 
where a is the average number of demands per unit time. 
The lead time, t, has a gamma type probability density 


function 


=ut ( ut ) k=1 


; _ ue 
f(tia,k) = ZIT 


k k 
with mean a and variance “ee 


Then the probability of exactly z demands during the lead time is 


+8 ] 
t) ue (ut) 
(zsk u a) -[ eat)” ° SU A at 
ae z! (k=1)! 
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letting x = (atu) t, then 
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23k,U,8) = pgp x Te eke 
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: atu . 
the last step being valid if k i: an integer. Thus the demand during 
lead time has a negative binomial distribution with 
mean = == 
and 


ak a 
variance = = (1 + 4 
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6, Summary and Conclusions 

The distribution of demand during lead time has been discussed 
when discrete random variables describe the demand and lead time. 
It is unfortunate that the distributions do not appear as well 
known distributions, The use of probability generating functions 
has been shown as a method by which the distribution can be found 
and used. Calculations of the probabilities are cumbersome since 


in the examples shown each of coefficients of the 3* 


in the prob= 
ability generating function involved an infinite sum, The amount 

of computer time required for the calculations has not been investi- 
gated, but since the calculations are straight forward and the number 
of terms required before the infinite sum converges to an approximte 
value are relatively small, the use of the procedure outlined in 
Section 4 should be useful and not too cumbersome to preclude their 
use. 

The use of the generating function method requires data for the 
estimation of the parameters in the demand and lead time distribu-— 
tions. This data is readily available in the historical files of 
an inventory system, Other methods of estimating the demand during 
lead time requires collection of data directly from the number of 
demands during a lead time. This data is not usually collected and 
even if it were available sufficient data to obtain a good statis- 
tical sample would require a very long time. 

Further investigation of distributions calculated by the methods 
in Section 4 may show that the probability distributions could be fit 
to some known distributions and the resulting distribution used to 


approximate the demand during lead time. 


Oe 





This investigation has been made primarily to consider a method 
of determining the demand during lead time distributions for low 
demand items. Low demand items are particularly troublesome in an 
inventory system, Military spare parts inventories have large 
numbers of very low demand items, The low demand items are aiso 
characterized by high variances [6]. 

The use of the Poisson and negative binomial distributions 
have been used extensively in this paper since these distributions 
have been found to fit low and medium demand items. The Poisson 
distribution is generally used for items with a demand of two or 
less items per year. The negative binomial distribution has been 
used for items with demands ranging from two to one hundred per 
year [7], In the case where the Poisson demand distribution is 
applicable, the generating function method for establishing a re= 
order point can be readily applied. Since the reorder point will be 
relatively low, the number of probabilities that must be calculated 
before the sum is equal to or greater than the service level desired 
will be small. It has not been discussed in this paper, but when 
reorder points have been previously established, the probability 
generating function for demand during lead time can be used to find 
the tail of the distribution. With a reorder point previously estab—- 
lished, a good estimate of the value of reorder point is available 
and with several trials on either side of the old reorder point, 
the tail of the distribution could be used to determine the new order 
point. In this manner, the probabilities would not have to be summed 


until the sum was greater than the desired probability. Feller [2] 


Be 





gives the following relations between the probability generating 
function and the generating function of the tail of the distri- 
bution. 


Let Q (s) = the generating function of the tail 


Then 


Q (s) _ Iz hls) 


les 


The calculation of the tail in this paper has not been investigated 
since it gets more complex than the probability generating function. 
However, in special cases and ina particular situation as described 
above, the tail of the distribution may be the easiest method for 
setting a reorder point. 

The special case of Poisson type demand and gamma distributed 
lead times was discussed because this combination of demand and lead 
time distributions lead to a negative binomial distribution of de- 
mand during lead time, Further investigation of this model with 
actual data may be very fruitful. The Poisson distribution is a 
good distribution for low demand items and the gamma distribution 
is a very versatile two parameter distribution that can be used to 
approximate a wide range of lead time patterns [5]. Consequently, 
any demand and lead time process that may possibly be fit to this 
model should be considered, especially since the probability law of 
demand during lead time may be readily computed, 

The writer wishes to thank Professor Robert R. Read of the 
U.S. Naval Postgraduate School for his suggestions, assistance, en= 


couragement and cooperation in the preparation of this paper. 
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APPENDIX I 
MEAN AND VARIANCE FROM PROBABILITY GENERATING FUNCTIONS 
The mean and.variance of a discrete random variable may be 
calculated from the probability generating function, It may be 


shown that (see e.g., Feller [ 2)) 
ete cate 1 ) 
Var(X) = f''(1) + £'(1) - fers}? 


As an example, the mean and variance of demand during lead 


time when lead time is a translated Poisson random variable 


a 





p(xsA) = en 


(x-1)! ee 


with mean = (A +1) and variance = A 


and demand has the negative binomial distribution 


p(x;r,p) = ( x q » X20,1,2,°°° 


with mean = a and variance = 4 
Pp Pp 


then the random variable Z has a probability generating function 


(ay ~ pie 


taking the first derivative with respect to s we have 


h'(s) = rqp* (e.\ exp Slee “yh + rape? ) (taf eefa2y 
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and 


h'(1) = sen Sy rq A 
pp 
xr 
=—- (14d) 
p 


the second derivative evaluated at one is 


h' t(4) = als } [rats 1° + rena + Ge } 
Pp 


then 


Var(z) = n'(1) + nt(1) - [ae(a}]? 


(=3-)"a4ay? + (Fe) 2. (414A) + a" (14A) = (2)? (140)" 


if 


= (44d) + [9 
p (; 


These results can be verified against the general formulas for 


the mean and variance given in Section 2. 
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